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This work presents a comparative analysis of hybrid systems based on molten carbonate fuel 
cells and making use of closed-cycle externally-heated bottoming systems. Two options are 
considered: reciprocating (Stirling) engines and supercritical carbon dioxide turbines. These 
two engines share the common feature of working on closed cycles with pure fluids (H 2 and 
C0 2 respectively) but, at the same time, they differ in their internal structure: Stirling engines 
make use of volumetric machinery whereas the supercritical carbon dioxide gas turbine 
(SC0 2 ) is composed by turbomachinery. In both cases, the working fluid is subjected to very 
high pressures and temperatures in the range of 50—200 bar and 40—650 °C. 

A brief description of both bottoming systems is provided in the first sections of this 
article along with their expected performance in on-design and off-design (part-load) 
operation. The analysis of each system is therefore split into three stages. First, the most 
relevant features of the models of performance are discussed. Then, a comparison is 
shown for on-design operation aiming to evaluate the maximum efficiency attainable by 
the proposed engines. Finally, an analysis of off-design operation is presented. 

The final section of paper concludes that hybrid systems based on atmospheric fuel 
cells and externally heated closed-cycle bottoming engines have the potential to outper¬ 
form conventional pressurised fuel cell plus gas turbine hybrids while avoiding the 
demanding operating conditions on the fuel cell topping system of the latter configuration. 
Copyright © 2012, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights 

reserved. 


1. Introduction 

The aim of this work is to develop a new concept of hybrid 
system based on topping molten carbonate fuel cells and 
bottoming externally-heated engines. The integration is of the 
indirect type (i.e. there is heat exchange but no mass transfer) 


which brings about a longer life expectancy of the fuel cell due 
to its smooth atmospheric operation in steady conditions. 
Ease of operation (load control) and the possibility to run the 
fuel cell in stand-alone mode are other interesting features. 

The first heat engine considered is a closed-cycle gas turbine 
operating with supercritical carbon dioxide (SC0 2 ). This is 
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a configuration studied extensively so far by the authors, who 
report a better performance of such system with respect to 
conventional hybrids (with either indirect or direct integration) in 
on-design and also off-design conditions [1,2]. The advantage of 
supercritical fluids is the very low work required to elevate the 
fluid’s pressure at the compressor and hence the higher useful 
work of the system. The interest of carbon dioxide lies in the close 
to ambient critical temperature (30.98 °C) and low critical pressure 
(73.8 bar), which make it possible to develop supercritical cycles 
with conventional cooling techniques and economical piping. 

The second heat engine under consideration is a Stirling 
engine. It is also a closed-cycle engine but makes use of 
a reciprocating configuration. The candidate working fluids of 
a Stirling engine are hydrogen, helium, nitrogen or simply air; 
in the present paper, hydrogen is selected for it yields the best 
performance [3]. 

A Stirling engine comprises: (i) a pair of cylinders where the 
fluid is compressed at low temperature and expanded at high 
temperature, (ii) a heater and a cooler to add and reject heat to 
and from the engine respectively, (iii) a regenerator located in 
the middle of these heat exchangers in order to increase the 
efficiency. The main advantages of this engine are compact¬ 
ness, low vibrations, low noise and good part-load perfor¬ 
mance. Durability and efficiency are also to be expected to 
depend on the design operating conditions [4], 

The benefits of incorporating either of these two engines 
into a high temperature fuel cell are discussed in this paper in 
a structured manner. First, conceptual overviews of the refer¬ 
ence fuel cell and hybrid systems are presented. Then, the two 
bottoming systems under consideration are described focusing 
on the modelling assumptions employed to estimate their 
performances. Finally, the efficiency in on-design and off- 
design operation is assessed, showing an important gain 
when supercritical C0 2 turbines are used instead of the 
conventional hot air turbines. 


2. Basic layout of the reference hybrid system 

A general layout of the reference hybrid system is shown in Fig. 1 
where a Molten Carbonate Fuel Cell is fed with preheated air and 
a mixture of water steam and partly-reformed natural gas. This 
raw fuel is reformed internally in order to take advantage of the 
heat released by the oxidation of hydrogen. The steam neces¬ 
sary for this process is generated outside the cell in a dedicated 
heat recovery steam generator that recuperates low 


temperature waste heat from the system. This option is 
preferred to recirculating a fraction of the cell exhaust gases 
because it allows reducing the stack temperature (hence 
increasing global efficiency) by optimising the temperature 
levels at which heat is exchanged by different subsystems of the 
plant. In other words, high temperature heat is used to drive the 
bottoming system whereas steam is generated at lower 
temperature. 

A catalytic burner is located downstream the cell exhaust to 
burn the inevitable excess fuel and increase the temperature of 
the heat source used to operate the bottoming engine (let it be 
noted that a minimum temperature of 700 °C is required if the 
heat engine is to achieve high efficiency). These hot gases enter 
a heat exchanger in which heat is transferred to the bottoming 
system (HX4 in Fig. 1) and, then, a fraction of them is recircu¬ 
lated to provide the cathode of the cell with the carbon dioxide 
necessary to avoid carbonate starvation in the electrolyte. The 
remaining gas is used to preheat the fuel and air streams into 
the cell and to generate steam for the internal reformer. 


3. Models of performance 

3.1. Fuel cell 


This works makes use of a lumped-volume model of perfor¬ 
mance of molten carbonate fuel cells with indirect internal 
reforming, of which a very detailed discussion is available in 
reference [5] by the authors. The fundamental equations that 
govern the operation of the fuel cell are listed below: 


• Nernst voltage equation: ideal open circuit potential 
depending on the pressure, temperature and composition of 
the reactants. Compositions are averaged logarithmically 
between the inlet and outlet sections of the cell. 


AG f _ AG° R U T n /WPoA 

2 -F 2 -F + 2 -F ^ P h 2 o J 


Ru-T 

2-F 
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/ Pc0 2; an 
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• Ohmic resistance: internal resistance to the flow of elec¬ 
tronic/ionic charge within electrodes/electrolyte. 


R ir = A ir -exp(AH ir /(R u -T)) (2) 

• Activation losses at the anode: responsible for the voltage 
drop due to the need to overcome the activation potential of 
the oxidation half-reactions. 


FUEL 



Fig. 1 - Hybrid system layout. 


Ran = A an -p H °' 5 -exp(AH an /(R u -T)) (3) 

• Activation and concentration losses at the cathode: the first 
term is conceptually similar to Eq. (3) whereas the second 
takes into account the rate of consumption of carbon dioxide 
to form carbonate ions at this electrode (C0 2 +0 2 +2e“ -► CO 3 ). 

Rea — ^ca,i"Po 2 •Pco 2 -exp(AH ca ,i/(R u -T)) 

+ -A ca2 ' Xco 2 ' exp(AH ca , 2 /(R u -T)) (4) 

• Operating current/voltage of the cell: Faraday’s law is used 
to evaluate the rate of consumption of hydrogen. This 
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information is then combined with the definition of fuel and 
carbon utilisations (U f and U C o 2 ) to determine the (raw) fuel 
flow rate and the rate of exhaust gas recirculation to the 
cathode inlet respectively. For the fuel, a standard hydro¬ 
carbon fuel C a Hp is considered, yielding the definition of the 
steam-to-carbon ratio given in Eq. (9). This parameter 
determines the amount of steam that has to be supplied to 
the reforming process to promote the conversion of hydro¬ 
carbons to hydrogen rather than the formation of elemen¬ 
tary carbon. It takes values in the range from 2 (the risk of 
carbonisation increases) to 5 (voltage drops due to the 
dilution of reactants at the anode). 


Table 1 - Compressibility factor Z for air and supercritical 
C0 2 at the inlet/outlet to/from the compressor/turbine. 

Location 

Variable 

Air 

s-co 2 

Compressor inlet 

Tr/Pr 

2.259/0.027 

1.014/1.017 


Z 

1 

0.25 

Compressor outlet 

Tr/Pr 

3.294/0.080 

1.271/3.050 


Z 

1 

0.60 

Turbine inlet 

Tr/Pr 

6.994/0.077 

3.037/2.929 


Z 

1 

1.02 

Turbine outlet 

Tr/Pr 

5.588/0.027 

2.609/1.059 


Z 

1 

1 


V = E 


Nernst 


— (Rtotal) J 






j ,cathode inlet 



Pi = (p| n - P° ut ) /In (p‘ n /p° ut ) (8) 

STCR = t - . nH2 ° . - (9) 

ftco + ftC0 2 + E • OL 


3.2. Supercritical carbon dioxide turbine 

The supercritical carbon dioxide power cycle has been 
extensively reviewed in literature in the last few years, start¬ 
ing from Dostal’s review [ 6 ] of the works by Angelino and 
Feher [7,8]. This growing interest has given way to dedicated 
technical meetings (Supercritical C0 2 Power Cycle Sympo¬ 
sium held once every two years in US) and a number of 
research works, most of which focus on nuclear power 
generation in high temperature gas reactors. More recently, 
SC0 2 systems have been proposed for solar applications [9] 
and hybrid systems with fuel cells [ 10 ]. 

Developing a model of performance of closed-cycle gas 
turbines involves the analysis of two types of components: 
turbomachinery (compressor and turbine) and heat 
exchangers (recuperator, heater and condenser/cooler). So far, 
different authors have approached the simulation of this 
turbomachinery based on the well founded thermodynamic 
and aerodynamic principles that are commonly used to design 
air compressors and gas or steam turbines [11,12]. Neverthe¬ 
less, this is not correct when applied to supercritical C0 2 
systems as noted by some works in the subject [13,14]. The 
markedly non-ideal behaviour of carbon dioxide near the 
critical point brings about the need to develop a specific 
database on which an efficient design of the compressor (and 
to a lesser extent the turbine) can rely. This is better assessed 
in Table 1 where the compressibility factors at the most 
relevant stations of the cycle are summarised for air and 
supercritical carbon dioxide . 4 

4 The compressibility factor Z is an index to evaluate how close 
to ideal the behaviour of a particular gas at a certain pressure and 
temperature is: p-u = Z-R-T. Thus, Z = 1 for ideal gases. 


This scenario suggests that a generalised performance 
map be employed in lieu of the analytical (equation-based) 
models typically used to estimate the performance of 
compressors in gas turbine engines. In order to have as 
accurate as possible estimations of the compressor perfor¬ 
mance, this map is adapted from the available experimental 
information gathered at the supercritical carbon dioxide test 
rig at Sandia National Laboratories, where a centrifugal 
compressor running on C0 2 at supercritical inlet conditions 
was tested in 2009 and 2010 [15]. For the turbine, the vanes are 
assumed to be choked and hence a single equivalent-nozzle 
curve is matched to the compressor map to yield the 
running line of the system in a pressure ratio vs corrected 
mass flow plot. The outcome of this approach is illustrated in 
Figs. 2 and 3 where the corrected mass flow rate in the 
compressor and turbine is defined as follows: 

m . . 

m c =- ( 10 ) 

Ro,in ' fl 0,in 


m t = 


fi^/f 0,in 


m 


O.in 



Note that subscript 0) in in Eqs. (10) and (11) indicates stag¬ 
nation conditions at the compressor/turbine inlet. It is also 
worth noting that the corrected mass flow is expressed 
differently in compressor and turbine, due to the different 
behaviour of the fluid at the inlet to each component. In effect, 



Fig. 2 - SC0 2 compressor performance map [15]. 
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Fig. 3 - SC0 2 turbine performance map [16]. 


Table 1 shows that carbon dioxide behaves close to ideally at 
the turbine (compressibility factor close to unity) whereas 
strong real gas effects are experienced at the compressor 
(compressibility factor ranging from 0.25 to 0.60). This latter 
feature of the fluid is accounted for by introducing R, y and p in 
the definition of m c as deduced from the application of Buck¬ 
ingham’s IT theorem [16]. 

Heat exchangers are modelled in this work by means of an 
e-NTU approach where effectiveness z and number of transfer 
units NTU are defined as follows: 

C* = C m in/C max = (™' C p) m i n /(hi • c p)max — ^ C^) 


^ Q. Chot ('fhot,in fhot,out) C co id (T co ld ou t T co id i n ) ( 13 ) 

Qmax C m in (Thot,in — 'f'colcLin) Lmin (Thot,in — T co id j n ) 


NTU = A-U/C min (14) 

C* in Eq. (12) is the heat capacity ratio where the heat 
capacity C of a fluid is defined as the product of mass flow 
times the specific heat at constant pressure. T stands for 
temperature of the cold (c) and hot (h) fluids at the inlet (in) 
and outlet (out) sections of the heat exchangers. Based on Eqs. 
(12) and (13), a specific correlation between z and NTU exists 
depending on the type and internal flow arrangement of the 
equipment. For a cross flow PFHE or PCHE, this equation is: 

= 1 - exp( - (1 - C*)NTU) 

1-C*exp(-(1-C*)NTU) 1 } 

The combination of Eqs. (14) and (15) provides an estimation 
of the global heat transfer coefficient and exchange area of the 
component. A more detailed description of the calculation 
procedure is available in reference [1]. 

Off-design operation of the gas turbine is performed by 
a combined system based on inventory and by-pass controls. 
At high loads, the output of the turbine is controlled by 
modifying the circulating mass flow in the system whilst the 
turbine inlet temperature remains constant. This is accom¬ 
plished by charging and discharging an auxiliary storage 
vessel connected in parallel to the compressor (Fig. 4). 
Nevertheless, this control strategy brings about an accelera¬ 
tion of the shaft at reduced loads which finds a limit at 
a certain rotating speed for which the centrifugal forces acting 
on the blade tips become too high. For further load reductions, 



Fig. 4 - Gas turbine layout showing control subsystems. 


a secondary control strategy is required which, in this case, is 
a turbine by-pass valve. This secondary control strategy 
consists on reducing the expanding mass-flow at the turbine 
while keeping it constant at the compressor. The net shaft 
work is so reduced. 

It has to be noted here that other control strategies like 
modifying the turbine inlet temperature are also possible and 
even yield fairly high efficiencies as shown in Fig. 5. Further¬ 
more, the comparison shown in this latter figure confirms 
that, apart from inventory control, the cycle is more efficient 
when a temperature control is adopted as opposed to using 
a turbine by-pass valve. Nevertheless, in spite of this obser¬ 
vation, using an inventory plus turbine inlet temperature 
combined control system is not feasible in practice since both 
strategies result in higher shaft speeds at reduced loads. That 
is to say, inventory control is limited by the highest shaft 
speed that the impeller/rotor tolerates. If a turbine inlet 
temperature were adopted for further load reduction, the 
shaft would speed-up bringing about a mechanical failure of 
the rotating parts of the C0 2 system. It must be remarked that 
this is because of the very high sensitivity of compressor 



Load SC02 f%l 


Fig. 5 - Part load performance of the SC0 2 cycle. 
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performance to inlet density (p 0 i) which varies dramatically 
when the low pressure level (p 0 i) of the cycle is reduced at 
part-load. Moreover, Eq. (13) shows that a density reduction 
would inevitably lead to higher corrected flow rates and, 
consequently, higher rotating speeds. Eventually, the 
maximum shaft speed for safe operation would be exceeded 
and the turbine would fail. Hence, an inventory followed by 
turbine by-pass system is most convenient. 

In light of the results above, the following control system is 
recommended as shown in Fig. 5: 

• 100—50% load range: inventory control. The lower limit of 
50% is set as the load for which the shaft rotates at 110% its 
rated speed. 

• Loads lower than 50%: by-pass control. The shaft speed 
reduces to about 80% of its rated value at 15% output. 


The reference engine is an alpha-type squared-four unit 
with double-acting cylinders, in which a 90° angle delay 
between each pair of compression-expansion volumes makes 
it possible to develop four working cycles in each crankshaft 
turn. One such pair of compression/expansion cylinders is 
shown in Fig. 7 where the internal geometry of the engine is 
clearly illustrated: compression cylinder (C), cooler (K), 
regenerator (R), heater (H) and expansion volume (E); inter¬ 
connection volumes between adjacent components are also 
accounted for in the calculations. 

The model is based on solving the equations of conserva¬ 
tion of mass, momentum and energy at each volume (node) 
shown in Fig. 7: 

• Mass conservation: 


3.3. 


Stirling engine 


dm/dt = rh m - rh in 


(17) 


The working cycle of a Stirling engine shares common 
features with that of the SC0 2 turbine since both are based on 
compressing a gas when it is cold and expanding it when it is 
hot. Nonetheless, the reciprocating layout of Stirling engines 
incorporates differential aspects as shown in Fig. 6 where the 
working cycle of these systems is depicted. The area enclosed 
by the p—V diagram represents the indicated (thermody¬ 
namic) work per cycle developed by a pair of cylinders: 


Wind = /p dv 


(16) 


The models of performance of Stirling engines are usually 
classified in categories depending on their complexity: zeroth 
order (experimental curve fitting), first order (rough estima¬ 
tions based on simplified models), second order (local effects 
like pressure drop or heat transfer are incorporated), third 
order (conservation of mass, momentum and energy applied 
to elements of the fluid domain within the engine) and fourth 
order (computational fluid dynamics applied). In this paper, 
a third order model of performance is used whose complete 
details are given in Ref. [17]. 





Energy conservation: 


= Q + rh in -(h in + 0.5-ufJ - m out -(h Pout + 0.5-u* ut ) - W 


(18) 


Momentum conservation: 


u-du + 


f-v 


2d 


dx + — = 0 


h 


(19) 


The previous equations are complemented by the following 
equation of state for real gases [18]. The utilisation of this 
equation makes it possible to account for the dense gas effects 
of hydrogen at high pressure and low temperature as follows 
(k = 1.9 x 10 -6 for hydrogen): 


P r V i = m r R u -(T i + fe-P 1 ) 


( 20 ) 


Starting from the principle of mass conservation applied to 
the engine as a whole, Eq. (21), and incorporating the equation 
of state just shown, Eq. (20), a codable equation to calculate 



Fig. 6 — Principle of operation of a Stirling engine. 


Fig. 7 — Calculation domain for the reference engine (pair of 
cylinders). Note local refinement of heat exchangers. 
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the reference pressure in the compression cylinder at time t is 
obtained, Eq. (22): 

n k Hr 

m c +m CK +m Kj + m KR +^m R . + m EH ++ m HE +m E = M 

i=l i=l i=l 

( 21 ) 


now studied. To this aim, a reference 500 kWe fuel cell is 
considered with the operating conditions summarised in 
Table 2. 

The design conditions of the supercritical carbon dioxide 
turbine are shown in Table 3. The most relevant information 
concerns the thermodynamic state of the fluid at the inlet to 



where F* is the ratio from the pressure in volume i to the 
reference pressure: 

Fi = Pi/P c (23) 

An iterative process is implemented in order to calculate the 
pressure inside the compression and expansion cylinders at 
each time step (and therefore the free volume inside the 
engine); the P—V diagram (and therefore the indicated work) is 
so obtained. The following assumptions are incorporated in 
the model: 

• For the calculation of pressure losses: 

o This effect is coupled to the thermodynamic model (i.e. it 
is not independent as it is in second order models). 

o Inertial losses and oscillatory flow effects are accounted 
for in the regenerator. On the contrary, the latter effects 
are neglected in circular pipes (heater/cooler) due to the 
high “fluid displacement/pipe diameter” ratio; stationary 
flow correlations are used instead [19]. 

o Pressure losses are also considered in the interconnection 
volumes. 

• For the heat transfer calculations: 

o Conductive heat transfer is accounted for within the gas 
(stream-wise) and between adjacent elements of the 
cooler, regenerator and heater walls. 

o Oscillatory flow correlations are used to evaluate Nusselt 
number in the regenerator whereas stationary flow 
correlations are used in circular pipes. 

The model so developed has been validated against 
experimental data available at the NASA Lewis Research 
Centre [20] showing a more than satisfactory agreement [17]. 


4. Hybrid system performance 

4.1. On-design operation 

The performances of hybrid systems comprising a topping 
molten carbonate fuel cell and either a bottoming supercrit¬ 
ical carbon dioxide turbine or a bottoming Stirling engine are 


the compressor, 308 K and 75 bar, which is very close to the 
critical point. For the specifications of turbomachinery and 
heat exchangers, these have been taken from an exhaustive 
literature review as reported by the authors in Ref. [2]. 

Finally, the rated operating conditions of the Stirling 
engine are listed in Table 4 where novel concepts to indicate 
the performance of the engine are needed due to its recipro¬ 
cating nature . 5 These parameters are the mean effective 
pressure (mep) and the mean pressure for mechanical losses 
(mpml). The mean effective pressure results from dividing the 
work developed (by a pair of cylinders) in a cycle by the 
volume displaced in the same period of time . 6 Similarly, the 
mean pressure for mechanical losses results from dividing the 
mechanical losses (experienced by a pair of cylinders) by the 
volume displaced in a cycle. These parameters are commonly 
used in the field of reciprocating engines and their mathe¬ 
matical descriptions are given below [ 21 ]: 

mep = f/P dv) j V D (24) 

mpml = Wioss/V D (25) 

A comparison of Tables 3 and 4 yields an interesting 
conclusion which is the higher efficiency of the supercritical 
carbon dioxide turbine when the same maximum tempera¬ 
tures are considered; this difference is in the order of 8—9%. 
Nevertheless, in spite of the lower efficiency of the Stirling 
engine, it is also remarkable that this latter engine is far more 
efficient than the conventional hot air turbine still proposed in 
many hybrid systems using atmospheric MCFCs. More in 
particular, the expected efficiency of such a system is around 
27%, as shown in Ref. [2] by the authors, what means that 
making use of a Stirling engine would provide a 30% gain in 
bottoming system efficiency. 


5 Note that the crankshaft mechanism impelling the flow of 
hydrogen from one cylinder to another makes it impossible to 
define a rated continuous flow like in a fuel cell or gas turbine 
engine. 

6 The mean effective pressure is a virtual pressure which acting 
on the piston from the top dead centre to the bottom dead centre 
(a complete stroke) develops the same effective work as the 
engine does in a cycle. 
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Table 2 - Fuel cell rated performance. 

Parameter 

Value 

Current density [A m 2 ] 

1100 

Voltage [V] 

0.716 

Anode/cathode pressure loss [%] 

1/2 

Fuel flow [g s -1 ] 

19.9 

STCR [-] 

3 

Temperature [K] 

923 

Fuel/C0 2 utilization [%] 

77.4/70 

Gross power output [kWe] 

500 

Auxiliary power [kWe] 

52.4 

Net efficiency [%] 

47.6 


These aspects are illustrated in Table 5 which provides 
a comparison between the rated performance of hybrid 
systems using either hot air turbine (standard case in the 
public domain), supercritical C0 2 turbines or Stirling engines. 
The following issues are to be noted: 

• The efficiency gain when a supercritical C0 2 turbine is used 
is remarkable. In particular, for the same waste heat avail¬ 
able from the topping system (MCFC), a 50% increase in 
output from the bottoming system is attained with respect 
to the same configuration using a hot air turbine. In terms of 
global power production, the previous effect has an impact 
of more than 5%. 

• The share of power delivered by the bottoming system 
when the SC0 2 turbine is used approaches 20%. 

• Using a bottoming Stirling engine in lieu of a hot air turbine 
yields a performance lying between that of a hot air turbine 
and a supercritical C0 2 system, even if it is closer to the first 
case. 

Another interesting aspect to be noted in Table 5 is the 
peak temperature of the Stirling engine, which is quoted as 
910 K. Given the volumetric nature of the system, the 
maximum temperature of the working gas within the engine 
is not constant, as opposed to the other turbine-based cases. 
Actually, this temperature depends on the volumetric flow of 
gases through the heater and on whether these gases flow 
from the compression cylinder to the expansion volume or 
vice versa. As a result, the peak temperature achieved by the 
working gas varies non-linearly in the range from 780 to 935 K, 


Table 3 - Rated performance of the SC0 2 cycle. 

Parameter 

Value 

Regenerator pressure loss hot/cold side [%] 

1.5/0.5 

Heater pressure loss [%] 

2 

Precooler pressure loss [%] 

1 

Compressor inlet pressure/temperature [bar/K] 

75/308 

Pressure ratio [-] 

3:1 

Turbine inlet temperature [K] 

923 

Compressor/turbine efficiency [%] 

90/95 

Recuperator effectiveness [%] 

95 

Mechanical efficiency [%] 

98 

Shaft work [kj kg -1 ] 

107.7 

Gross efficiency [%] 

39.1 


as a function of the crank-angle. Thus, an average value of 
910 K is indicated in Table 5. 

The performance of the system in part-load operation is 
illustrated in Fig. 8 which confirms that the utilisation of 
externally-heated closed-cycle bottoming systems is benefi¬ 
cial in terms of system efficiency. In particular, an efficiency 
gain of around six percentages is to be expected for a 50% load 
reduction in both cases (58-65% for the SC0 2 systems and 56 
to 61 for the Stirling-based hybrid). The root cause for this 
behaviour is found in the superb part-load performance of the 
bottoming cycles whose efficiencies remain fairly constant 
even at reduced load settings. This is particularly evident in 
the case of the supercritical C0 2 turbine, which exhibits a drop 
in efficiency from roughly 40% to around 37% in the afore¬ 
mentioned range. 

The information provided in Table 5 and Fig. 8 confirms that 
the integration of molten carbonate fuel cells and closed-cycle 
heat engines has the potential to boost the efficiency of 
atmospheric hybrid systems to values in the order of 60% 
while, at the same time, preserving the safe operating condi¬ 
tions of the cell. This is thanks to the variety of part-load 
mechanisms available for closed cycles, which yield 
a comparative advantage over conventional systems where 
temperature and shaft-speed control only are possible. Even 
though such control strategies demand for more complex 
subsystems (for instance inventory tanks or additional valves), 
the increased power output or reduced fuel consumption and 
the already complex layout of an internal reforming fuel cell 
system justifies the utilisation of the proposed configuration. 

Finally, some additional remarks must be highlighted with 
regard to the choice of supercritical carbon dioxide turbines or 
Stirling engines. First, despite the higher efficiency of the SC0 2 
option, Stirling engines are in a more advanced state of 
development though it must be acknowledge that they still 
lack economic competitiveness due to the very high 
manufacturing and maintenance costs. In this regard, it 
would be possible to have pilot MCFC-Stirling pilot plants 
installed in the short term. In contrast, SC0 2 turbines are still 
an immature technology as deduced from the very few 
experimental tests in operation which, in most cases, 
concentrate on one aspect of the system or on demonstrating 
its technical feasibility rather than on developing a product to 


Table 4 - Rated performance of the Stirling engine. 3 


Parameter Value 


Mean pressure [bar] 

150 

Shaft speed [rpm] 

1750 

Heater wall temperature [K] 

993 

Cooling water inlet temperature [K] 

323 

Effective work [kj cycle -1 ] 

976.5 

Mean effective pressure [bar] 

18.6 

Mean pressure for mechanical losses [bar] 

1.1 

Gross efficiency [%] 

36.07 


a Note that even though the effective work per cycle is given for 
the reference engine used for validation (30-35 kWe rated power 
output), the mean effective pressure has the noteworthy advantage 
of not being affected by the engine displacement. The same applies 
to the mean pressure for mechanical losses. 
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Table 5 - Rated performance of the hybrid systems. 

Parameter 

Air 

sco 2 

Stirling 

Effective power MCFC [kW] 

447.6 



Effective power heat engine [kW] 

71.7 

103.5 

79.8 

Global effective power [kW] 

519.2 

551.1 

526.7 

Ratio Power MCFC/Engine [%] 

86.2/13.8 

81.2/18.8 

85.0/15.0 

Efficiency MCFC [%] 

47.6 



Engine peak temperature [K] 

948.6 

947.1 

910 

Net efficiency bottoming engine [%] 

27.6 

39.8 

30.3 

Global net efficiency [%] 

55.2 

58.7 

56.0 



Hybrid system load [%] 


65 

60 

55 

50 
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40 

35 

30 

25 

20 

15 


v° 

>s 

o 

c 

a; 

o 

it= 

0) 

v 

o 

>s 

O 


Fig. 8 — Part-load efficiency of the hybrid systems under 
consideration. 


be deployed to the market. It is hence expected that these 
systems approach the commercial stage in the mid-term. 

Further to these considerations, another differential aspect 
between both bottoming systems is the output range. Hence, 
Stirling engines are commercially available in the range of 
tens of kilowatts (<100 kW), with some companies reporting 
to have plans to develop units able to deliver hundreds of 
kilowatts (200—500 kW) [22], whereas the supercritical carbon 
dioxide turbine is expected to be cost-effective for power 
outputs approaching 1 MW and beyond [23]. In light of these 
considerations, Stirling engines are most convenient for 
small-scale hybrid systems whilst SC0 2 turbines should be 
adopted in mid and large-scale hybrid power plants. 


5. Conclusions 

This paper presents a comprehensive assessment of the most 
interesting bottoming system to be used in hybrid systems 
composed by high temperature molten carbonate fuel cells 
and bottoming heat engines of the closed cycle type. The 
rationale behind this analysis is to explore the benefits that 
could eventually be derived from the flexible features of 
closed power cycles as opposed to open systems: working 
fluid, system pressure, part-load strategy and others. Two 
candidate engines are considered: Stirling engines 


(volumetric) and closed-cycle turbines operating on super¬ 
critical carbon dioxide. 

The main conclusions drawn from this work are: 

• Comprehensive and accurate modelling tools are provided 
for each subsystem (fuel cell, Stirling engines and SC0 2 
turbine), with the capacity to simulate the performance in 
on-design and off-design operation of each one of them. 

• The utilisation of a supercritical carbon dioxide turbine is 
expected to boost the efficiency of the corresponding hybrid 
system to values in the order of 60% for systems comprising 
fuel cells operating at atmospheric pressure. 

• Using a Stirling engine is somewhat less interesting due to 
the limited efficiency gain with respect to standard hybrid 
systems that incorporate hot air turbines. The main reason 
for this is the effect of the volumetric nature of the unit on 
the peak temperature achieved by the working gas. This 
effect brings about a reduction in engine efficiency. 

• The efficiency of the SC0 2 turbine remains almost unaltered in 
a wide range of load-settings, hence providing a remarkable 
part-load performance of the corresponding hybrid system. 

• This pattern is somehow observed in the Stirling-based 
configuration though the performance of this system is worse. 

• SC0 2 turbines and Stirling engines fit to different output ranges. 
Hence, Stirling engines are suitable for outputs lower than 
a hundred kilowatts whilst SC0 2 turbines are most interesting 
for outputs in the range of a few hundred kilowatts or more. 

Further research work by the authors is aimed at assessing 
the costs of both systems and the potential to benefit from 
economies of scale, even if such an analysis incorporates 
important uncertainty. Once this complementary information 
is available, the technical comparison provided in this paper 
would be integrated with the economic information in order 
to identify the most cost-effective technology to be used in the 
next generation of hybrid systems. 

REFERENCES 


[1] Sanchez D, Chacartegui R, Munoz de Escalona JM, Sanchez T. 
Performance analysis of a MCFC & supercritical carbon 
dioxide hybrid system under part load operation. Int J 
Hydrogen Energy 2011;36:10327-36. 

[2] Sanchez D, Chacartegui R, Munoz de Escalona JM, Sanchez T. 
A comparison between molten carbonate fuel cells based 
hybrid systems using air and supercritical carbon dioxide 
Brayton cycles with state of the art technology. J Power 
Sources 2011;196:4347-54. 

[3] Karabulut H, £inar C, Oztiirk E, Yiicesu HS. Torque and 
power characteristics of a helium charged Stirling engines 
with a lever controlled displacer mechanism. Renew Energy 
2010;35:138-43. 

[4] Kongtragool B, Wongwises S. A review of solar-powered 
Stirling engines and low temperature differential Stirling 
engines. Renew Sust Energy Rev 2003;7:131-54. 

[5] Munoz de Escalona JM, Sanchez D, Chacartegui R, Sanchez T. 
A step-by-step methodology to construct a model of 
performance of molten carbonate fuel cells with internal 
reforming. Int J Hydrogen Energy 2011;36:15739-51. 

[6] Dostal V, Driscoll M, Hejzlar P. A supercritical carbon dioxide 
cycle for next generation nuclear reactors. Report No. 












578 


INTERNATIONAL JOURNAL OF HYDROGEN ENERGY 38 (2013) 57O-578 


MITANPTR-100. Cambridge: Massachusetts Institute of 
Technology; 2004. 

[7] Angelino G. Carbon dioxide condensation cycles for power 
production. Paper 68-GT-23. ASME; 1968. 

[8] Feher EG. The supercritical thermodynamic power cycle. In: 
Proceedings of the intersociety energy conversion 
engineering conference; 1967. Miami Beach. 

[9] Ma Z, Turchi C. S-C02 power cycle configurations for 
concentrating solar power systems. In: Proceedings of the 
2nd CSP plant optimization conference and expo; 2011. San 
Francisco. 

[10] Sanchez D, Chacartegui R, Jimenez-Espadafor F, Sanchez T. A 
new concept for high temperature fuel cell hybrid systems 
using supercritical carbon dioxide. J Fuel Cell Sci Technol 
2009;6:02136. 

[11] Gong Y, Carstens NA, Driscoll MJ, Matthews IA. Analysis of 
radial compressors options for supercritical C02 power 
conversion cycles. Topical Report MIT-GFR-034. Cambridge: 
MIT; 2006. 

[12] Moisseytsev A, Sienicki JJ. Development of a plant dynamics 
computer code for analysis of a supercritical carbon dioxide 
Brayton cycle energy converter coupled to a natural 
circulation lead-cooled fast reactor. Report ANL-06/27. 
Argonne National Laboratory; 2006. 

[13] Wang Y, Guenette G, Hejzlar P, Driscoll M. Compressor 
design for the supercritical C02 Brayton cycle (Paper 
AIAA2004-5722). In: Proceedings of the 2nd IECEC 
conference; 2004. Rhode Island. 

[14] Ulizar I, Pilidis P. A semiclosed-cycle gas turbine with carbon 
dioxide-argon as working fluid. J Eng Gas Turb Pow 1997;119: 
612-6. 

[15] Wright S, Radel RF, Vernon ME, Rochau GE, Pickard PS. 
Operation and analysis of a supercritical C02 Brayton cycle. 
Report SAND2010-0171; 2010. Alburquerque. 

[16] Sanchez T, Munoz A, Jimenez-Espadafor FJ. Turbomaquinas 
termicas. Madrid: Sintesis; 2004 (in Spanish). 

[17] Munoz de Escalona JM, Sanchez D, Chacartegui R, Sanchez T. 
Model of performance of Stirling engines. In: Proceedings of 
the ASME turbo expo; 2012. Copenhagen. 

[18] Vargaftix NB. Tables on the thermophysical properties of 
liquids and gases. 2nd ed. John Wiley a Sons, Inc.; 1975. 

[19] Tew RC, Thieme LG, Dudenhoefer JE. Recent Stirling engine 
loss-understanding results. Report NASA TM-103122; 1994. 

[20] Lundholm SGK, Moodysson BO. The design and development of 
a 40 kW Stirling engine-the power system. IMech E 1982;C27/82. 

[21] Heywood JB. Internal combustion engine fundamentals. 
Singapore: McGraw-Hill; 1998. 

[22] Corporate information Artemis Energy, www.artemis-energy.eu. 

[23] Sienicki J, Moisseytsev A, Fuller B. Scale dependencies of 
supercritical C02 Brayton cycle technologies and the optimal 
size for a next-step supercritical C02 cycle demonstration. 
In: Proceedings of the 2nd SC02 power cycle symposium, 
BIECEC conference; 2011. Boulder, CO. 


Glossary 

Latin 

a: sound speed, m s -1 
A: area, m 2 

C: heat capacity, W K” 1 


C p : specific heat, J kg” 1 K” 1 
d h : hydraulic diameter, m 
E: Nernst voltage, V 

/: friction factor, — 

F: Faraday constant, 96,485 C 
G: Gibbs free energy, J kg” 1 
h: enthalpy, J kg” 1 
j: current density, A m 2 
K: equilibrium constant, — 
m: mass, kg 

mep: mean effective pressure, bar 

mpml: mean pressure for mechanical losses, bar 

m: mass flow, kg s” 1 

M: total mass, kg 

MCFC: molten carbonate fuel cell 

NTU: number of transfer units, — 

n: molar flow, mol s” 1 

p: partial pressure, bar 

P: pressure, bar 

Q: heat, W 

R u : universal gas constant, 8.3145 J mol -1 K” 1 
R: resistance, Q m 2 

SC0 2 : supercritical carbon dioxide turbine 
STCR: steam to carbon ratio, — 
t: time, s 

T: temperature, K 
u: internal energy, J kg” 1 

U: overall heat transfer coefficient, W m~ 2 K” 1 
U C o 2 ■ carbon utilisation factor, % 

Uf: fuel utilization, — 
u: velocity, m s” 1 
V: voltage, V 
Vi.- volume, m 3 
W; work, J kg 1 
W: power, W 
x: molar fraction, — 

Z: polarisation, V 

Greek 

e: effectiveness, % 
p: density, kg m” 3 

Subscripts 
an: anode 

C: compression cylinder 
ca: cathode 

CK: connecting duct compression cylinder—cooler 
D: displaced 
E: expansion cylinder 
H: heater 

HE: connecting duct heater—expansion cylinder 
in :inlet 
ir: Ohmic 
K: cooler 

KR: connecting duct cooler—regenerator 

max: maximum 

min: minimum 

out: outlet 

R: regenerator 

RH: connecting duct regenartor—heater 



